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Abstract: We investigated the phytosociology, structure and dynamics 
of Pinus roxburghii in 40 stands in northern areas of Pakistan by using 
cluster analysis (Ward’s agglomerative clustering) and ordination (Non¬ 
metric Multidimensional Scaling). Cluster analysis revealed three major 
groups associated with specific environmental characteristics: (1) P. 
roxburghii (2) Pinus-Quercus baloot and (3) Pinus-Olea ferruginea 
community types. NMS-ordination showed the major gradient as an 
amalgam of elevation (r 2 = 0.441 ,p <0.01) and slope ( r 2 = 0.391, p <0.05) 
as the two topographic factors correlated with species distribution. The 
first ordination axis also showed positive correlation with soil variables 
like pH and electrical conductivity, suggesting that soil chemistry was 
related to topographic characteristics and probably acted as a secondary 
gradient. We also examined size class distributions, age structures and 
growth rates of the three communities in order to describe community 
development and dynamics. Total tree density was 14700 plants/ha, with 
P. roxburghii having a relative density of 82% to 100%. Density of juve¬ 
nile and total density and basal area of the subordinate tree species were 
low. The low density of trees in the smallest diameter size-class sug¬ 
gested that the recruitment of small P. roxburghii plants into the adult 
population may be lower than the required replacement rate for the 
stands. Pooled size-class distributions for the species showed a multimo¬ 
dal pattern with some regeneration gaps. Browsing, heavy logging and 
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other anthropogenic activities were the overriding factors responsible for 
the poor recruitment of P. roxburghii. We concluded from the age struc¬ 
ture that the forests were characterized by the dominance of young trees. 
Growth rate analysis revealed that P. roxburghii was the fastest growing 
species among the conifers species in Pakistan. In view of its relatively 
fast growth and longevity, P. roxburghii seems to be a suitable choice for 
short-term cultural practices in order to enhance wood production in 
lesser Himalaya and Hindukush ranges of Pakistan. 

Keywords: phytosociology, structure, dynamics, multivariate techniques, 
dendrochronology, Himalayan range 

Introduction 

In Pakistan, about 4.8% of the total area is forested which is very 
low in comparison to 23-70% in other countries of the world 
(Ahmed et al. 2010; Khan 2012). One-third of the forested area 
in Pakistan is covered by productive forests and the remaining 
two-thirds are categorized as protected forests (Sethi 2001). The 
sub-tropical pine forests located in the lesser Himalayan and 
Hindukush range of Pakistan are dominated by broad-leaved and 
coniferous tree species have significant contribution to produc¬ 
tive forest types. To date, eight conifer tree species have been 
reported (Nasir and Ali 1972) of which two species i.e. P. brutia 
and P. halepensis are not native to Pakistan. Among the native 
tree species Pinus roxburghii (Chir pine) is considered to be an 
important species (Gupta and Dass 2007) occurring at elevations 
from 900 to 1800 m above sea level (Siddiqui et al. 2009). In 
Pakistan this species is distributed in the lower parts of Khyber 
Pakhtunkhwa, Azad Jammu and Kashmir in sub-tropical dry and 
moist temperate areas (Ahmed et al. 2006), while it is also 
planted in certain areas of Punjab along roadsides and gardens. It 
is the only pine in the country of the northern forests with a small 
overlap with Pinus wallichiana (Kail) at the upper elevation limit 
(Saddozai 1995). P. roxburghii has a wide ecological amplitude 
and considerable economic importance, providing large stretches 
of grazing lands due to its typically well-developed grass layer 
(Wahab 2011) and valuable timber-wood and resin. Hence most 
of the population in its zone of occurrence depends on these 
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forests for various purposes (Gupta and Dass 2007). The eco¬ 
logical role of P. roxburghii forests is also important in terns of 
watershed, microclimate, anti-erosion, sanitary-hygienic, recrea¬ 
tion and conservation functions (Malabadi and Nataraja 2006; 
Ravindra et al. 2007). 

In addition to its ethno-botanical importance, P. roxburghii is 
also used as charcoal, dye, herbicide, ink, lighting, resin and 
wood, while the importance of these forests as a source of many 
other bio-resources including mushrooms, medicinally important 
herbs, shrubs and game can also not be ignored (Ghildiyal et al. 
2010). Owing to the economic importance of P. roxburghii, its 
intensive usage and direct devastation have gradually led to 
shrinkage of the forested area and degradation of this ecologi¬ 
cally and commercially important species (Nafeesa et al. 2007; 
Miehe et al. 2009). Most of these stands are isolated and have 
reduced wood properties, low productivity and decreased farm¬ 
ing and resource significance. The large scale destruction of 
forests in many sectors has resulted in several environmental 
problems (Singh 1985; Ives and Messerli 1989) that have caused 
changes in vegetation structure, diversity and composition (Gaur 
1982). Knowledge of the current floristic composition and struc¬ 
tural attributes is useful in understanding the prevailing status of 
commercially and ecologically important species because this 
underpins resource description, land-use planning and conserva¬ 
tion management (Singh 1976; Kumar et al. 2004). 

The structural and functional parameters of P. roxburghii as¬ 
sociations in Pakistan were reported by Ahmed et al. (2006). 
Information related to these communities was also published by 
other workers (e.g. Ashraf 1995; Ajaib and Khan 2004; Ajaib et 
al. 2004; Malik et al. 2007; Nafeesa et al. 2007; Siddiqui et al. 
2009; Shaheen et al. 2011; Wahab 2011). Some dendroecological 
and dendrochronological studies have been reported by Ahmed 
et al. (2009). However, little is known of the phytosociology, 
structure and dynamics of the species. Therefore, the present 
study was designed with the following specific objectives: (1) 
quantify the relationship of species composition and stand struc¬ 
ture attributes of P. roxburghii forests with relation to major 
environmental variables in northern Pakistan, (2) elucidate 
growth and development patterns of the forest species with rela¬ 
tion to site variables, (3) compare the pine forest in the study 
area with other forests on national and regional levels. Results 
described here may be used as reference conditions for similarly 
situated forests on the lesser Himalaya and Hindukush range that 
are being managed for old growth characteristics ( i.e. growth rate, 
age, diameter and vertical structure). 

Materials and methods 

Study area and climate 

The northern areas of Pakistan cover about 72500 km 2 spanning 
from 34° to 37° N and 71° to 74° E (Ahmed et al. 2011). There 
are three well-known ranges of mountains i.e. the great Himala¬ 
yas (average elevation 4000 m), which stretch for about 2500 km 
from east to west and border the Hindukush range (5000 m) to 
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the northwest and Karakorum to the north (Sethi 2001). Geo¬ 
graphically, the area is linked with China through the famous 
Silk Road which is roughly 4700 m above sea level and reported 
to be the highest sealed highway in the world (Brian 1991). The 
area adjoins the disputed territory of Kashmir to the east and 
Afghanistan to the west (Sheikh 1993). The lesser Himalaya is of 
middle elevation (1800-4500 m), represented by the Pir Panjal 
range, and increases in elevation from south to north. The Hin¬ 
dukush range has a number of passes i.e. Shandur pass (3700 m) 
that connect Gilgit-Baltistan and Chitral (Geoffreyl984), 
Shangla pass (3000 m) links Swat valley to the upper parts of the 
Indus valley while the Lowari pass (3118 m) connects Chitral to 
Dir, Swat and Peshawar (Manfred et al. 1986; 2000). 

Major portions of these mountains are generally bare of vege¬ 
tation though rich forests are found in some areas (Champion et 
al. 1965). However, P. roxburghii is completely absent in the 
Karakorum range (personal observations). These mountain rang¬ 
es have a variety of micro-climates typically of highland type 
(Fowler and Archer 2006) and consequently support a variety of 
vegetation types (Shaheen et al. 2011). Winters are long, cold 
and snowy while summers are short and mild. Rainfall is highest 
in highland zones and usually associated with elevation (Hussain 
et al. 2010), while some areas like north-west Gilgit and Chitral 
experience low rainfall as they are rain shadow areas beyond the 
reach of the monsoon (Khan et al. 2013). 

Field methods 

Vegetation survey and soil analysis 

After the general reconnaissance of the study area, 40 stands 
dominated by P. roxburghii were sampled throughout their natu¬ 
ral limits by Point Centered Quarter (PCQ) method (Cottam and 
Curtis 1956). The sampled forest stands were located both in 
Hindukush and Himalayan ranges of Pakistan, 2-3 ha in total 
area, and showed no sign of recent major disturbances. Understo¬ 
rey vegetation including herbs, shrubs, seedlings and saplings 
was also sampled in 3 m x 3 m quadrats following Ahmed and 
Shaukat (2012). Twenty-five PCQ points and quadrats were 
sampled at regular intervals of 20 m in each stand as suggested 
by Ogden and Powell (1979). Prior to collection of the vegeta¬ 
tion data, information was recorded including slope angle (Sunto 
clinometer), aspect (compass), elevation (altimeter), geographi¬ 
cal coordinates (GPS) and the presence of disturbances or human 
interference. Diameter of all trees (DBH > 10 cm; height 1.3 m 
above ground) and height of saplings (DBH < 5 cm) were meas¬ 
ured following Ogden and Ahmed (1989). Seedlings (height less 
than 30 cm) were counted and all recorded species were pre¬ 
served and identified by reference to literature in Pakistan (Nasir 
and Ali 1972). Soil samples were collected to a depth of 25 cm 
using a soil auger at three different locations within a stand and 
samples were pooled to obtain a composite soil sample. The 
composite soil samples (CSS) from each stand were subjected to 
physicochemical analysis following the procedures described by 
Sparks (1996) and Carter and Gregorich (2008). Soils were 
sieved through a 2-mm sieve in order to separate pebbles and 
gravel. Water holding capacity was determined in accordance 
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with the measurements of Keen (1931) and Carter and Gregorich 
(2008). Soil organic matter (%) was quantified by loss on igni¬ 
tion at 500°C. Soil pH was measured in a 1:5 soil: distilled water 
paste. Major cations (Ca i+ , Mg 2+ , K + , Na + ) were quantified by an 
atomic absorption spectrophotometer (Pie-Unicam) after extract¬ 
ing in 0.5 M ammonium acetate. Total soil nitrogen was esti¬ 
mated by the Kjeldahl method (Carter and Gregorich 2008). 

Tree-ring width data 

Dendrochronology was studied following the general methodol¬ 
ogy described by Fritts (1976) and Speer (2010). The tree sam¬ 
ples in the fonn of cores were taken using a Swedish increment 
borer from healthy and sound trees that were also free from se¬ 
vere competition in accordance with the practice of Cook and 
Kairiukstis (1990). Attempt was made to extract two core radii 
from each living tree at breast height above the ground, parallel 
to the slope contour as suggested by Fritts and Wang (1986). The 
extracted samples were kept in plastic straws with both ends 
sealed for safety and were labeled with the relevant information 
(Huang and Zhang 2007). The core samples were dried at room 
temperature and mounted in a wood skeleton using wood glue 
following Heinrich (2004). The surfaces of the core samples 
were smoothed using progressively fine grades of sandpaper 
following Bowers (1964) and Pilcher (1990). 

Laboratory methods 

Data analysis 

Importance value and absolute values (density/ha and basal area 
m 2 /ha) were computed following Curtis & McIntosh (1950) and 
Mueller-Dombois and Ellenberg (1974). Densities of herbs, 
shrubs, seedlings and saplings were calculated on per ha basis. 
For effective analysis of the vegetation and related environ¬ 
mental variables, both classification and ordination techniques 
were employed. The input data matrix of tree vegetation of 40 
forest stands based on importance value was developed and sub¬ 
jected to cluster analysis and ordination using PC-ORD version 
5.10 (McCune and Grace 2002). We applied Ward’s agglomera- 
tive clustering technique (Orloci and Kenkel 1985) for numerical 
classification of tree species and the results were presented in a 
dendrogram. No species was excluded as an outlier from the 
analysis due to the low number of species. 

Non-metric multidimensional scaling (NMS) was chosen as 
the ordination method to investigate the underlying trends in 
vegetation distribution patterns of the tree species and to corre¬ 
late the compositional patterns with the set of environmental 
variables. NMS ordination is often preferred over other popular 
methods because of its non-linear non-parametric basis, which 
provides effective ordination of vegetation data where species 
distributions usually follow non-linear (often Gaussian distribu¬ 
tion) patterns (Shaukat 1989; Enright et al. 2005). This ordina¬ 
tion method avoids the assumption of linear relationships among 
variables and the use of rank distances tends to linearize the rela¬ 
tionship between distances measured in environmental space. 
This relieves the “zero-truncation problem” that plagues all ordi¬ 
nations of heterogeneous community data sets (Minchin 1987; 


McCune and Grace 2002). It allows the use of any distance 
measure or relativization (MuCune and Mefford 2005). Enright 
et al. (2005) suggested the use and selection of 1, 2 and 3 dimen¬ 
sions for best solution of interpretation, in terns of trade-off 
between complexity associated with increased dimensionality 
and reduction of stress function. We selected the best solution for 
each dimensionality by choosing the option of the lowest stress 
value (Timilsina et al. 2007) along with the initial configuration 
of correspondence analysis (Shaukat 1989). The relationships 
between environmental variables and NMS ordination of stands 
were examined by superimposing and contouring the stand data 
on the ordination configuration. Variables used for NMS were 
topographic, soil and vegetation variables. The topographic vari¬ 
ables were elevation and slope, while soil variables included 
organic matter (OM %), pH, Carbon (%), Nitrogen (%) Magne¬ 
sium (%), Sodium (%) Potassium (%) and total dissolved salts. 
Statistical analyses (ANOVA, Pearson product moment correla¬ 
tion and regression analysis) were performed using SPSS Ver. 14 
and Sigma Plot while the group mean vectors of environmental 
variables were tested by the multivariate statistic “Hoteling’s T 2 ” 
using the program HOTELLING developed in C~ + . 

Size class distribution 

The diameter of P. roxburghii was used to construct size-class 
distributions and associated tree species were not included. 
Measurements for P. roxburghii trees from all 40 stands were 
pooled following Ahmed et al. (2011) into three groups because 
of low density. The pooled data were divided into a series of 
classes, i.e. 10-19.9 cm, 20-29.9 cm. For multi-stemmed trees, a 
diameter was assigned that was equivalent to a hypothetical sin¬ 
gle stem with a cross-sectional area equal to the sum of all the 
cross-sectional areas of the individual stems (Russell and Fowler 
1999; Ryniker et al. 2006). Juveniles of P. roxburghii less than 
150 cm in height or less than 3 cm in circumference were ex¬ 
cluded from size class structure. 

Tree-ring data processing 

All cores were not included in the age and growth rate analysis 
due to missing portions caused by stem-rot in the center of some 
trees: only those cores crossing the center of the trunk were used 
in the analysis. Individual cores were examined under a variable 
power binocular microscope for visual cross-dating (Strokes and 
Smiley 1968). The age and growth rates of the core samples were 
calculated following Roza (2004) and Lv and Zhang (2012). In 
cases where cores did not pass through the center or pith of the 
tree, we followed Ogden (1980) to calculate the number of miss¬ 
ing rings and add these to the total age of the tree. Three cross- 
sections were taken at ground level of saplings in each stand and 
the rings on these sections were counted. It was assumed that the 
average height and growth rate shown by these saplings could be 
used to approximate the time required for the tree to reach the 
height at which cores could be taken (Ogden 1981; Ahmed et al. 
2009). These years were added to the age of each core to obtain 
the total age of the tree. The length of the cores was divided by 
the number of rings present in the core and the average growth 
rate in years per centimeter was calculated. Linear regression 

Springer 




514 


Journal of Forestry Research (2014) 25(3): 511-521 


was computed for each site to establish the relationships between 
age, growth rate and tree diameter. For better age estimation, the 
average age of two cores of the same tree were employed. 


Results 


Forest structure 


Ward’s agglomerative cluster analysis resulted in three major 
groups with 75% of the information in the abundance of species 
retained (Fig. 1). The percent chaining in the cluster analysis was 
10.44%, indicating the effectiveness of cluster analysis. Three 
major groups/associations were identified as described below. 


25 


50 


75 


100 


Group-I 


Group-II 


di 




St.l 
St .2 
St.3 
St.5 
St .6 
St.7 
St.8 
St.ll 
St.12 
St. 13 
St.14 
St. 15 
St. 16 
St. 17 
St.18 
St .20 
St.21 
St.22 
St .2 3 
St.24 
St.25 
St.27 
St.28 
St.29 
St.30 
St.31 
St.32 
St .4 
St. 10 
St.19 
St.9 
St.26 
St.34 
St.3 6 
St.3 7 
St.3 5 
St.33 
St.38 
St.40 
St.39 


Fig. 1: Dendrogram showing the different associations identified by the 
hierarchical agglomerative cluster analysis based on importance value of 
trees. 

P. roxburghii association 

Group I, the P. roxburghii association, was represented by 27 
sites/stands (Fig. 1). Among these sites P. roxburghii occurred as 
a pure community (100% IVI) and the trees were distributed 
exclusively at dry and moist sites of the subtropical areas of the 
Hindukush and Himalayan ranges (Table 1). The range of den¬ 
sity ha" 1 in this pure association was 41—1490 trees/ha with a 
mean of 398±69 trees/ha (Table 2). The total density in this 
group was 14700 trees/ha. The mean density of seedlings 
(234±51 individuals/ha) was higher than that of saplings (166±34 
individuals/ha) in this group. 


Table 1: Mean importance value in the three groups of tree vegetation 
derived from Ward's cluster analysis. 


Tree species _ 

Group 1 

Group 2 

Group 3 


Mean±SE 

Mean±SE 

Mean±SE 

P. roxburghii 

100±00 

62±2.8 

62.142.6 

Q. baloot 

X 

32±4.9 

0.31440.2 

Q. dilitata 

X 

5±4.91 

2.8441.0 

A. modesta 

X 

X 

4.141.3 

B. papyriffera 

X 

X 

1.740.43 

C. serrate 

X 

X 

0.540.6 

M. buxifolia 

X 

X 

8.544.9 

O. ferruginea 

X 

X 

17.543.6 

P. granatum 

X 

X 

241.0 

F. palmata 

X 

X 

0.540.3 

Table 2: Mean density/ha in the three groups derived from Ward’s clus- 

ter analysis. 




Tree species _ 

Group 1 

Group 2 

Group 3 


Mean±SE 

Mean±SE 

Mean±SE 

P. roxburghii 

398±69 

85±9.25 

386478.99 

Q. baloot 

X 

24±10.96 

1342.5 

Q. dilitata 

X 

16±5.89 

3.541.0 

A. modesta 

X 

X 

1643.9 

B. papyriffera 

X 

X 

8400 

C. serrate 

X 

X 

2400 

M. buxifolia 

X 

X 

28412.5 

O. ferruginea 

X 

X 

53415.7 

P. granatum 

X 

X 

541.8 

F. palmata 

X 

X 

1400 

Table 3: Mean basal area m 2 /ha of the three 

groups derived from Ward’s 

cluster analysis. 




Tree species _ 

Group 1 

Group 2 

Group 3 


Mean±SE 

Mean±SE 

Mean±SE 

P. roxburghii 

129±33 

57±16.09 

125432 

Q. baloot 

X 

10±4.50 

441.6 

Q. dilitata 

X 

6±1.99 

1.440.8 

A. modesta 

X 

X 

3.140.9 

B. papyriffera 

X 

X 

0.2400 

C. serrate 

X 

X 

0.02400 

M. buxifolia 

X 

X 

6.40.85 

O. ferruginea 

X 

X 

1442.88 

P. granatum 

X 

X 

240.37 

F. palmata 

X 

X 

0.32400 

Note: x is absence of species 

in different groups. 


Basal area of the dominant species ranged from 9- 

-898.9 nr/ha 

with a mean of 129±33 nr/ha (Table 3), while the total value for 

basal area was 3468.99 

m 2 /ha. A number of herbs and shrubs 

were recorded apart from the seedlings and saplings of the over- 

story tree species. Among these Dodonea viscosa 

, Indigofera 
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gerardiana, Buddleja crispa, Isodon rugosus, Maytenus 
royleanus , Heteropogon contortus, Viola biflora and Micromeria 
biflora were frequent species with frequency ranging from 
10-100% in this group. 

P. roxburghii-Quercus baloot association 

Group II, the P. roxburghii-Q. baloot association, was recorded 
at three sites. Two sites in this association were from the sub¬ 
tropical moist temperate area and one site represented the dry 
temperate area. Among the tree species P. roxburghii was domi¬ 
nant with 62±2.8% importance value followed by Q. baloot 
(IV=32±4.9%) and Quercus dilatata (IV=5±4.91%) (Table 1). 
Mean tree density of the leading species was 85±9.25 trees/ha 
(range 66-98 trees/ha), whereas the total trees numbered 253 
trees/ha, less than for groups I and III (Table 2). The basal area 
ranged from 38 to 89 m 2 /ha with a mean of 57±16.09 m 2 /ha. The 
overall basal area of P. roxburghii was 171 m 2 /ha. The density of 

Q. baloot ranged from 3 to 40 trees/ha (of a total of 72 trees/ha 
of all species) with an average of 24±10.9 trees/ha. Q. baloot 
basal area ranged from 0.87 to 15.5 m 2 /ha (of the total for all 
species of 29.23 m 2 /ha) with a mean value of 10±4.50 m 2 /ha. Q. 
dilatata was a sub-ordinate species at a density of 16±5.89 
trees/ha with an average basal area of 6±1.99 m 2 /ha (Table 3). 
Saplings of the dominant tree were more numerous as compared 
to the total saplings of the other species (12-33 individual ha" 1 ). 
Dodonea viscosa, Berberis lycium, Mallotus philippensis, Myrtus 
communis, Viola biflora and Fragaria nubicola were dominant 
understory herbs and shrubs at 20 to 70% frequency. 

P. roxburghii-Olea ferruginea association 
Group III, the P. roxburghii-Olea ferruginea association, was 
found at 10 sampled locations. These were located in the dry 
temperate area though a few sites were near moist temperate sites. 
Ten tree species were recorded in the group in which P. rox¬ 
burghii was highly dominant (IV=62.1±2.6%) followed by broad 
leaved O. ferruginea (IV=17.5±3.6%) and Monotheca buxifolia 
(IV=8.5±4.9%) (Table 1). Three species, Acacia modesta, Quer¬ 
cus dilatata and Punica granatum, were comparatively abundant, 
whereas the remaining species had importance values of less 
than 1%. The overall density and basal area of trees in this group 
was 2521 individuals/ha and 1400 m 2 /ha, respectively. P. rox¬ 
burghii shared a total density of 1251 trees/ha (aver- 
age=386±78.99 tree/ha) and a total basal area of 435.77 m 2 /ha 
(average=125±32 m 2 /ha). Among the broad leaved species, O. 
ferruginea had the highest density (average=53±15.7 trees/ha) 
and basal area (average=14±2.88 m 2 /ha) followed by M. buxifo¬ 
lia at an average density of 28±12.5 trees/ha and basal area of 
6±0.85 m 2 /ha. Other species had average densities less than 8 
trees/ha with a small fraction of the basal area (Tables 3 and 4). 
Group III, however, had highest species richness followed by 
Group I. In addition, density values for juveniles were higher for 
the dominant species followed by Olea ferruginea and M. buxifo¬ 
lia. The understory was dominated by D. viscosa at 100% fre¬ 
quency followed by Teucrium stocksianum (50%), Ajuga par- 
viflora (50%), Indigofera gerardiana (30%) and Otostegia lim- 
bata at 20% frequency. Phagnalon niveum, Filago hurdwarica. 


Limonium cabulicum and Periploca aphylla were the sparsely 
distributed species. 


Table 4: Means of environmental variables for the groups derived from 
Ward’s cluster analysis. 


Factors 

Group 1 

Group 2 

Group 3 

Mean±SE 

Mean±SE 

Mean±SE 

Elevation 

1322T43 

1137±75 

1734T481 

Slope 

31±2.1 

28±3.8 

19±1.0 

pH 

5.2±0.08 

6±0.44 

6.3±0.17 

WHC % 

56±1.2 

48±1.0 

53±1.16 

Salinity 

0.15±0.01 

0.30±0.08 

0.26T0.03 

Cond. 

230±22 

237±22 

307±65 

TDS 

167±17 

120±11 

265±34 

OM % 

8.3±0.39 

6.8T0.5 

7.7±0.43 

Ca % 

0.40±0.06 

0.37±0.15 

0.45±0.09 

Mg % 

0.45±0.02 

0.61±0.02 

0.54T0.04 

Na % 

0.77±0.09 

0.61±0.14 

0.55T0.03 

K % 

1.6±0.11 

1.9T0.23 

1.6T0.20 

N % 

0.19±0.06 

0.14±0.04 

0.17±0.03 


Note: WHC= water holding capacity %, OM= organic matter, TDS= total 
dissolved solids 

Environmental characteristics of the groups 

The environmental variables recorded for the three groups 
showed marked differentiation (Table 3). Group I was set apart 
with low elevation, low slope angle, low pH, and high water 
holding capacity (WHC %), organic matter (OM %), total dis¬ 
solved salts (TDS), sodium (Na %), magnesium (Mg +2 ), potas¬ 
sium (K %) and nitrogen (N %). Calcium (Ca +2 ) concentrations 
were higher than at Group II sites and lower than at Group III 
sites. Group II was characterized by low elevation, water holding 
capacity, organic matter, salinity, conductivity, and Nitrogen. 
Conversely, slope gradient, pH and Mg +2 content were higher in 
this group. Finally, Group III was characterized by high eleva¬ 
tion, low slope angle and high pH and Ca +2 content. The salient 
feature of this group also included high conductivity, salinity and 
K + content in comparison with the other groups. 

Statistical analysis 

The groups derived from Ward’s clustering technique with their 
corresponding environmental variables were analyzed by uni¬ 
variate analysis of variance (ANOVA). Elevation and slope were 
significantly correlated (F — 7.5; p <0.01). The edaphic variables 
including soil pH, water holding capacity, salinity and conductiv¬ 
ity were non-significant. Soil organic matter and nutrients (Ca, 
Mg, Na, K and total N) showed non-significant differences be¬ 
tween group means. 

The multivariate assessments of ecological parameters of the 
groups were also compared based on mean vectors by the appli¬ 
cation of Hotelling’s T 2 test. The topographic, edaphic and soil 
nutrients were compared between pairs of groups with respect to 
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the first set of variables (topographic). The comparison between 
group I and II yield a significant difference (F =3.4, df =27, df 2 
=10, p <0.05). However, the mean vector for groups I and III 
were not significant (F =1.43, df\ = 27, df 2 =3 ns), while groups 
II and III had significantly different mean vectors (F =5.8, dfi 
=10, df 2 =3 p <0.01). With respect to the edaphic variables, the 
comparison of groups I and II showed non-significant difference 
(F =0.9, dfi =27, df 2 =10 ns). Similarly, groups I and III were 
similar (F =0.462, dfi =27, df 2 =3 ns), as were groups II and III 
(F =0.487, dfi=10, df 2 =3 ns). Mean vectors for nutrients in 
groups I and II, I and III and II and III did not differ significantly. 

Forest-environment relationships (Ordination) 

The 2-D NMS ordination of tree vegetation data (stands ordina¬ 
tion) is shown in Fig. 2. The results of stand data yielded a two 
dimensional configuration with a comparatively high (20%) 
stress value. The groups/community types derived from Ward’s 
agglomerative clustering strategy were superimposed and were 
clearly separated in the ordination space. The results revealed 
that axis 1 alienated the P. roxburghii pure community from the 
P. roxburghii - Q. baloot and P. roxburghii - O. ferruginea com¬ 
munity types. P. roxburghii pure community type was distin¬ 
guished from the other two groups on axis 2. 



Fig. 2 : NMS ordination of tree vegetation (stands) representing three 
major groups derived from cluster analysis were superimposed on the 2- 
D ordination space. 

The correlation of environmental variables and the NMS ordi¬ 
nation axes suggested a strong relationship between elevation ( r 2 
=0.441, p <0.01) and slope (r 2 =0.391, p <0.05) along the first 
ordination axis. Soil pH was also significantly correlated (r 2 
=0.423, p <0.01) with ordination axis 1 (Table 5). Salinity, con¬ 
ductivity and total dissolved salts were also significantly corre¬ 
lated with ordination axes 1 and 2 (p <0.05). Among the other 
measured variables, water holding capacity (WHC), organic 
matter (OM) and available Ca ++ , +2 , Mg +2 , K + , Na~ and N + did 
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not yield significant relationships with any of the ordination axes 
(Table 5). The environmental parameters were also tested using 
one-way ANOVA and indicated statistically significant differ¬ 
ences between the groups in tenns of topographical variables. 


Table 5: Pearson’s product moment correlation coefficients between the 
ordination axes and the environmental variables examined. 


No. 

Variables 

Axis 1 

Prob. Level 

Axis 2 

Prob. Level 

1 

Elevation 

0.441** 

p<0.01 

0.181 

Ns 

2 

Slope 

0.391* 

p< 0.05 

0.191 

Ns 

3 

pH 

0.423** 

p<0.01 

0.042 

Ns 

4 

WHC % 

0.139 

Ns 

0.042 

Ns 

5 

Salinity 

0.315* 

p< 0.05 

0.334* 

p< 0.05 

6 

Cond. 

0.326* 

p< 0.05 

0.358* 

p< 0.05 

7 

TDS 

0.322* 

p< 0.05 

0.337* 

p< 0.05 

8 

OM % 

0.302* 

p< 0.05 

0.211 

Ns 

9 

Ca 

0.081 

Ns 

0.025 

Ns 

10 

Mg 

0.064 

Ns 

0.151 

Ns 

11 

Na 

0.057 

Ns 

0.118 

Ns 

12 

K 

0.09 

Ns 

0.056 

Ns 

13 

N 

0.083 

Ns 

0.134 

Ns 


Significance levels: *: p <0.05; **: p <0.01 and Ns: Non-significant 
Size class structure 

Size class structure based on diameter of trees across the 40 
stands was pooled into three groups (Fig. 3). The pooled diame¬ 
ter distribution showed an uneven size structure, with monotoni- 
cally decreasing density with increasing tree size up to the largest 
DBH class in group I, which comprised a secondary peak in 
density. The peak in the largest DBH class resulted in part from 
inclusion in the class of all trees greater than 40 cm. P. rox¬ 
burghii was present in all DBH classes but well represented in 
small and medium DBH size classes. M. buxifolia was absent 
from the lowest size class, but was present in other DBH classes 
at low numbers. The most abundant species in the 10-20 cm 
DBH class was O. ferruginea (47%) followed by A. modesta and 
Q. baloot. P. granatum were more abundant in < 30 cm DBH 
classes and rare in classes above 30 cm in all forest stands. These 
species exhibited an unstable population size structure. 

Size class distribution of group II was completely different 
than for groups I and III. Group II showed a typically uneven 
pattern of DBH size classes with some regeneration gaps within 
the larger size classes (Fig. 3b). In this group, P. roxburghii rep¬ 
resented only 33% of trees in the small size classes while in the 
medium size classes the density ha' 1 of trees was higher (50%) 
than in the small and large size classes. Generally, the number of 
trees was negligible and regeneration gaps were a prominent 
feature of this group, particularly in the larger size classes. The 
size frequency distribution in group III (Fig. 3c) revealed a pat¬ 
tern almost comparable with that of group I. In this group the 
tree ha' 1 in smaller and medium size classes was 47% and 26%, 
respectively. However, nearly 14% of P. roxburghii were in 
large and 12% in giant size classes. 
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r 2 =0.9436) while, age and growth rate and diameter and growth 
rate were also significantly correlated (Table 7). 


Table 6: DBH, age and growth rates of P. roxburghii in three groups 
derived from Ward's cluster analysis 


Groups 

Species 

Diameter 

(DBH cm) 

Age 

Growth rate 

(year/cm) 

Tree sam¬ 
ples 

i 

P roxburghii 

57±2.9 

115±6.5 

2.5±0.09 

N=125 


Range 

12-155 

20-360 

0.78-5.2 


2 

P.roxburghii 

4143.3 

54±5.7 

2.8±0.09 

N=60 


Range 

10-113 

15-198 

1.2-5.3 


3 

P roxburghii 

45±2.9 

69±5.4 

2.9±0.13 

N=74 


Range 

11-116 

14-209 

1.1-5.5 



Table 7: Regression analysis of DBH, age and growth rates of P. rox- 


burghii 

in three groups derived from cluster analysis. 



Groups 

Equations 

(DBH vs 

Age) 

R 2 

DBH vs 

Growth rate 

R 2 

Age vs 

Growth 

rate 

R 2 

i 

y= 20185x + 0.7597*** 

y= 

0.0589 Ns 

y= 

0.1269 


0.047 


0.0035x4 


0.012x4 


2 

y= 1.663 lx 

0.9436*** 

y= 0.0117x 

0.9436*** 

y= 0.197x 

0.4599* 


+ - 14.034 


4 2.256 


4 2.0836 


3 

y= 1.7458X 

0.9065*** 

y= 0.0176 x 

0.4945* 

y= 

0.4398* 


- 10.074 


4 1.7288 


0.0304x 4 



Significance levels: *: p <0.05; ***: p <0.001 and Ns: Non-significant 

P. roxburghii in group III had diameters ranging from 11-116 
cm (45±2.9, n =74) and age from 14-209 years (69±5.4, n =74). 
Regression analysis showed a significant positive relationship 
between diameter and age (y =1.7458x - 10.074, r = 0.906). 
Increment of growth ranged from 1.1-5.5 with an average of 
2.9±0.13 years/cm. Similar to group II, the trees in groups III 
exhibited significant linear relationship between age and diame¬ 
ter (r 2 =0.4945, p <0.005) and diameter and growth rate (r 2 
=0.4398 ,p <0.05) respectively (Table 7). 


Fig. 3a, b and c: Pooled size class structure of Pinus roxburghii based on 
groups obtained from Ward’s cluster analysis of 40 stands. Diameter 
classes are 10-19.9, 20-29.9, 30-39.9 cm and so on. 

Age and growth rates 

The diameter of P. roxburghii trees in group I ranged from 
12-155 cm (Mean ± SE= 57±2.9) while the age ranged from 
20-360 years (115±6.8, n=125). Age was positively correlated 
with diameter (r 2 = 0.7597, p < 0.001). Annual increment was 
estimated to range from 0.78-5.2 years/cm (2.5±0.09, n=125) 
and showed non-significant linear relationships with diameter 
and age (Table 6). Diameter of trees in group 2 ranged from 
10-113 cm (41±3.3, n = 60) and ages ranged from 15 to 198 
years with a mean annual increment of 2.8±0.09 years/cm. Tree 
diameter was significantly correlated with age (y =1.663x-14.034, 


Discussion 

The classification and ordination techniques showed well- 
defined patterns in vegetation composition and yielded comple¬ 
mentary results. Ward’s agglomerative clustering procedure 
partitioned the vegetation into three main groups based on the 
importance value of tree species. In our study this division was 
made on the basis of species composition for the entire study 
area. P. roxburghii was the dominant species in terns of its im¬ 
portance value in all stands because it typically grows in pure 
stands throughout Pakistan, India, Nepal and Bhutan (Sangye 
2005; Gupta and Dass 2007; Ghildiyal et al. 2009; Siddiqui et al. 
2009). The members of Fagaceae, Sapotaceae and Moraceae are 
often in association as dominant species on lower elevations 
(Khan et al. 2011). However, with P. roxburghii the members of 
these families co-occur rarely and whenever associated they are 
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represented with very low importance values. Species like Q. 
baloot, Q. incana and M. buxifolia are the common, almost ex¬ 
clusive low altitude evergreen representatives that are found in 
association with P. roxburghii (Badshah et al. 2010; Khan et al. 
2010). Singh and Singh (1986) reported that Q. incana and Q. 
baloot are widely distributed in the west with higher concentra¬ 
tion in the central Himalaya. However, in the present study both 
species were poorly distributed in few of the forest stands sur¬ 
veyed though frequently distributed at lower elevations in the 
Hindukush and Himalayan regions (Khan 2012). Q. dilatata is 
typically a west Himalayan species but is poorly represented in 
Pakistan (Ohsawa et al. 1986; Sheikh and Kumar 2009) and was 
absent from the entire study sites. Generally, Quercus and other 
broad leaved species occurred with low importance values in the 
sampled communities. 

NMS-ordination of arboreal vegetation (stand ordination) de¬ 
picted a vegetation continuum that appeared to be a function of 
altitude and slope as both are correlated with the first ordination 
axis. Thus topographic gradient was the predominant environ¬ 
mental gradient that controlled the composition of the vegetation 
as has been reported in other portions of the Hindukush and Hi¬ 
malayan mountains (e.g., Saddiqui et al. 2010; Ahmed et al. 
2011; Khan et al. 2013). The first ordination axis also showed 
positive correlation with soil variables like pH and electrical 
conductivity. Apparently, the variables related to soil chemistry 
were associated to topographic characteristics and they probably 
acted as a secondary gradient (cf. Greig-Smith 1983) which per¬ 
haps influences the uptake of nutrients. Such soil patterns appear 
general for mountain slopes where temperature decreases and 
rainfall increases with elevation (Allen and Pet 1990). 

The density-ha" 1 of P. roxburghii and associated broad leaves 
species in the present study corresponds with the low altitude 
species composition recorded for the lesser Himalayan and Hin¬ 
dukush range of Pakistan. The density values in the present study 
are within the range values that are reported by Ahmed et al. 
(2006), Wahab et al. (2008), Ahmed et al. (2009), Siddiqui et al. 
(2009) and Khan et al. (2010) in the low-land dry and moist 
broad leaved and evergreen conifer forests of the Hindukush and 
Himalayan ranges. The density of P. roxburghii varied signifi¬ 
cantly between the present study sites. Density was substantially 
lower at dry sites than at moist sites and P. roxburghii occurred 
more frequently on north facing slopes. Basal area (nr-ha" 1 ) re¬ 
corded in this study ranged from 9 to 898.9 nr-ha" 1 which is 
comparatively high for the dominant species and might be attrib¬ 
uted to high numbers of individuals in large size classes. How¬ 
ever, the basal area values for the associated broad leaved spe¬ 
cies were low in comparison to the range of 21.33-27.5 nr-ha" 1 
reported by other workers (Ahmed et al. 2009; Khan et al. 2010; 
Khan et al. 2011), from undisturbed forests in the lesser Hindu¬ 
kush. The associated broad leaved species contributed less than 
10% of the total basal area of all species. The low density and 
basal area for the associated broad leaved species were due to the 
low density of sizable trees, perhaps its harvesting for fuel-wood 
and because of the occurrence of P. roxburghii mostly in pure 
stands. 


The reverse J-shaped population structure of P. roxburghii in 
some stands in the study area is in conformity with many other 
forest stands in lesser Hindukush (Siddiqui et al. 2009; Wahab 
2011) and Himalayan forests (Ahmed and Naqvi 2005; Hussain 
et al. 2010;). The reverse J-shaped pattern for DBH frequency 
indicated good regeneration potential (Nath et al. 2005). This 
type of pattern explains indigenous factors and exogenic distur¬ 
bances (Barker and Kirkpatrick 1994). The negative exponential 
pattern is characteristic of undisturbed forests (Elouard et al. 
1997). However, in the current study P. roxburghii was the ma¬ 
jor contributor in lower size classes, hence showed maximum 
regeneration potential. The low representation of P. roxburghii in 
both sapling and seedling classes was observed, indicating com¬ 
paratively low level of regeneration. Generally, the presence of 
more small-sized tree than large-sized ones is considered as an 
indication of a regenerating forest (Bhat et al. 2011) but lower 
number of juveniles is treated as an indicator of the degenerating 
status of the stand (Condit et al. 1998). In the present study this 
phenomenon might be attributed to anthropogenic disturbances 
like illicit cutting and browsing of livestock. Absence of trees 
from particular size classes was also observed, indicating that 
mature trees were preferred for cutting or logging for beams and 
other purposes. Moreover, removal of trees and mortality create 
gaps and reduced over-crowding resulting quick growth of re¬ 
generating species that tend to be straight and tall (Obiri et al. 
2002). Such phenomena were observed in P. roxburghii stands 
that lead to the selective harvesting for specific uses from a par¬ 
ticular DBH class and non-selective tree cutting for fuel wood 
have been reported as affecting the forest composition and struc¬ 
ture (Saddiqui et al. 2009; Ahmed et al. 2011). The removal of 
the high number of trees in particular size class suggests that 
people need them in bulk quantity to meet their various needs. 
These practices in these forest sites are common due to its low 
elevation and amidst human habitation lead to rapid deterioration 
of these economically and ecologically important stands. 

The age and growth rate of the species suggested that this spe¬ 
cies is fast growing and relatively short-lived. In view of its fast 
growth and longevity, P. roxburghii seems a suitable choice for 
short-term cultural practices to enhance wood production in the 
lesser Himalaya and Hindukush ranges. These findings are also 
in agreement with Ahmed et al. (2009) and Ali (2011) who ob¬ 
served this species to be a favorite plant for re and afforestation 
purposes in the area. The examination of ring-width characteris¬ 
tics suggest that P. roxburghii may not be a suitable species in 
prediction of climatic changes, due to a number of reasons i.e. 
double ring/missing ring fonnation, excesses resin and fast 
growth. In this context, it is extremely difficult to cross-date, 
though studies have shown the dendroclimatic potential of the 
species growing at diverse ecological sites in western Himalaya 
(e.g. Pant and Borgaonkar 1984; Bhattacharyya et al. 1992; 
Schmidt 1993; Chaudhary et al. 1999; Cook et al. 2003). This 
difference may be due to the topography and other environ¬ 
mental factors that influence the growth of the species. We also 
concluded that diameter is a good predictor of age and growth 
rates in case of P. roxburghii as, both the parameters showed 
strong correlation with diameter. These findings accord well with 
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those of Ahmed and Sarangzai (1991) who reported that age 
increases with increasing diameter. It was reported that species 
from moist temperate areas generally grow faster (Ahmed et al. 
2009). However, such a statement is not applicable for all the 
species. In case of P. roxburghii from dry and moist areas, 
growth is faster than other conifer species because it prefers to 
grow on moderate to gentle slopes and flat areas mostly on the 
northern aspect. Therefore, faster growth may be dependent upon 
the better soil moisture regime of the stands rather than other 
extrinsic factors or its intrinsic characteristics. 

Based on these results, legal categorization and prohibition of 
extraction by the local communities has little effect on anthropo¬ 
genic pressure on forests, in deciding extraction pressure are 
village size, distance, accessibility and approachability. Unless a 
suitable management package is evolved and implemented, and 
conservation measures are taken to reduce and control the an¬ 
thropogenic activities, under the current regime, loss of forest 
degradation will be continued in northern Pakistan. 
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